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Abstract 

The treatment of a quantized field in a curved spacetime requires the study 
of backreaction of the field on the spacetime via the semiclassical Einstein 
equation. We consider a free scalar field in spatially flat Robertson- Walker 
space time. We require the state of the field to allow for a renormalized 
semiclassical stress tensor. We calculate the sigularities of the stress tensor 
restricted to equal times in agreement with the usual renormalization pre- 
scription for Hadamard states to perform an explicit renormalization. The 
dynamical system for the Robertson Walker scale parameter a(t) coupled to 
the scalar field is finally derived for the case of conformal and also general 
coupling. 



1 Introduction 

The studies of quantized fields in curved spacetimes usually assume a fixed back- 
ground spacetime on which a quantized field is defined [5j. Such a setting has 
also been used to investigate the simplest cosmological model of a homogeneous 
isotropic universe. In this context the mode spectrum of a scalar field in a homo- 
geneous isotropic universe undergoing an era of inflation has been found to be the 
famous and experimentally confirmed scale free Harrison-Zeldovich spectrum [9]. 

The work which is concerned with the backreaction of the field on the evolution 
of the spacetime itself usually takes a mean field approach in which quantum fluctu- 
ations only contribute to the effective potential of the (classical) expectation value 
of the field [10, 19J. The dynamics of the quantum degrees of freedom is not con- 
sidered. In other approaches, the quantum degrees of freedom either decouple from 
the spacetime [15] due to conformal coupling for a massless field, or a large mass 
of the field is assumed |3J, which implies another semiclassical approximation. This 
approximation in fact implies that the field configuration is dominated by renor- 
malization ambiguities. Again the coupling of quantum degrees of freedom to the 
geometry of the space time is only approximate. 
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In this work, which is partially based on we determine the coupling of a 
free quantum field to the scale parameter of the spatially flat Robert son- Walker 
spacetime. The result is an equation for a dynamical system with infinitely many 
degrees of freedom that can, at least in principle, be solved. 

The dynamical system is derived via an expansion of Riemann normal coordi- 
nates in the Robert son- Walker canonical coordinates up to third order [5]. This 
allows to calculate the singular terms of the Hadamard bidistribution and its first 
and second time derivative restricted to equal time surfaces up to constant terms 
(zero mode terms) of the second time derivative. Those terms are needed to renor- 
malize the energy-momentum tensor of the free field restricted to equal time sur- 
faces. The latter establishes the dynamics of the Robertson- Walker scale factor (up 
to renormalization ambiguities). 

While finalizing this paper we got aware of the publication [13], which derives 
exactly the same equation of motion as we do for the case of conformal coupling 
using a completely independent argument. In fact, this work in the conformally 
coupled case also establishes the existence of solutions for small time intervals given 
that the initial state fulfills certain conditions, which goes beyond the scope of this 
article. 



2 General description of dynamics and renormal- 
ization approach 

Throughout the article we restrict to flat Robertson- Walker (RW) spacetimes with 
space-time dimension 4. This allows to use the standard Fourier transform on spa- 
tial sections of constant RW time in order to formulate the dynamics of the field in 
terms of modes. As the stress tensor is formulated as a differential operator acting 
on the two-point function, we will formulate the dynamics of the field in such a way 
that the two-point function of the field and its associate momentum field operator 
restricted to the given time is the dynamical variable. We will also restrict our atten- 
tion to homogeneous, isotropic, quasifree and pure states that either are Hadamard 
states or are sufficiently close to Hadamard in the sense that they allow the same 
renormalization prescription for their energy-momentum tensor. Hadamard states 
and adiabatic vacuum states on RW space time have been studied e.g. by Liiders 
and Roberts Juncker and Schrohe [7] and more recently by Olbermann |12] . 

The line element of a spatially flat Robertson- Walker metric for a homogeneous, 
isotropic spacetime is 

d s 2 = dt 2 - a 2 {t)dx 2 (1) 

where we call a(t) the scale parameter and define the Hubble parameter H(t) = -. 
We consider a scalar free field of mass m coupled to the scalar curvature R{t) = 

-6 (H(t) + 2H 2 {t)^ with coupling £, 

(□-£i? + m 2 )0 = O, (2) 
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where □ is the D'Alembertian. 
the field modes (f>k is given by 



We use the abbreviation uj\ = a 2 k 2 + m? — £R. The Klein-Gordon equation for 



+ 3H<p k + u 2 k <p k = (3) 



>k 

Using the canonical momentum field Tik = a 3 <fik we consider the Hamiltonian form 
of this equation 

a -3 \ ( d> h 



It is shown in [UJ that the equal time two-point function, i.e. the two-point 
function on a Cauchy surface, of a state can be described by matrices 



(5) 



where = G^t, k) is the Fourier transform in z = x — y of two equal time 

field operators G(t, x, y) = uj((f)(t, x)<f)(t, y)), G^^ the Fourier transform of one equal 
time field operator and one canonically conjugated momentum operator etc., and k 
is the modulus of the momentum conjugated to z. Our normalization convention 
for the Fourier transform is JP(f)(k) = J R3 f(x)e lk ' x dx. The positivity of the state 
will enforce that the matrix is positive semidefinite with its determinant vanishing 
for a pure quasifree state. 

The symmetric part of the two-point function modes then fulfills the linear sys- 
tem of equations 



Gcf>(t>,k } 




I o 


2a~ 3 





G(07r),fc 


\-\ 


-<™* 





a- 3 


GftTttk J 




V o 


-2a 3 LU 2 k 







d t GW fe = -a 6 u 2 k a G (<fr),k (6) 



Here and in the following () stands for symmetrization and [] for anti sym- 
metrization in the field and momentum operator. This system of equations has one 
conserved quantity per mode 

Jk = G^^G^.k — G^ )k (7) 

which reduces the number of degrees of freedom per mode to two, as required. 

The condition for a state to be pure and to induce a representation of the canon- 
ical commutation (CCR) algebra then is 

\/k : J k = -Gf MJfc = \ (8) 

which implies the vanishing of the above-mentioned determinant for all modes. 
Next we turn to the semi classical Einstein equation. This equation has the form 

Gfiv = SnG (T^ (9) 
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where the expectation value of the stress tensor must be renormalized. In this 
equation is the Einstein tensor and G is the gravitational constant and both 
should not be confused with G^k etc. which are two-point functions. 

We will apply the point-splitting procedure as formulated in [T3], see also [16J. 
In order to preserve general covariance of the renormalization prescription, we need 
to subtract the Hadamard bidistribution fj, described in some detail in [S], from 
the two-point function before removing the point splitting. Actually, a sufficiently 
precise approximation S) n of the Hadamard parametrix does the job as well 

tfM) wA ,= I™ ^(x,y)[(0(x)0(y)) w -i3n(x,y)]+V(t;). (10) 
with , cf [EJUS], 

t^( v ) = S J U ^ J {Am 4 + Bm 2 R + CR 2 + DR a(3 R af3 ) d g x (11) 

where A, B, C and D are real valued renormalization degrees of freedom and d g x 
stands for the canonical volume form ^J\g\dx. 

D$. (x, y) is the symmetrization in x and y of the following second order partial 
differential operator (cf. [I4J Eq. (10) for the details, where we corrected some minor 
misprints) 

d'xA - h^(9 J '%A,s ~™ 2 )+i [{R, u - \g,»R) + 2g, v (□, + g^ff^) 

-2 {D' x ^ v + ff x M] - r]g fll/ (D x -£R+ m 2 ) (12) 

Here & = 8"(v,x)d XtV , S(v,x) being the geodesic transport and D' analogously 
denotes a covariant derivative. The quantities g, R and the Ricci tensor R pv are 

evaluated at v. Note here we use sign convention (H ) instead of ( — h ++) 

in jH] which flips the sign in the last expression in the brackets (and some more in 
the text below). 

The definition of S) n for d — 4 is 

1 u(x,y) 1 / a(x,y) 



4n 2 a(x,y) 4n 2 



log (-^jr) E y) k - ( 13 ) 

^ ' k=0 



Here a(x,y) is the squared geodesic distance, u(x,y) is the square root of the Van 
Vleck-Morette determinant (which for dimension 4 is Uo{x,y) as defined in |T3] up 
to normalization). u(x,y) and the Vk{x,y) (corresponding to Uk+i(x,y) of [14] ) are 
functions that can be determined by a recursive system of differential equations that 
depends exclusively on invariant quantities and the operator □ — £R + m 2 : 

2g> xv {x){d^a){d%u) + (D x a - A)u = (14) 
2g^(x)(d*a)(dZv ) + (D x a - 2)v = -(□* + m 2 - £R(x))u (15) 
2g^{x){d x a){d x u v k+1 ) + (D x a + 2k)v k+1 = -(D x + m 2 - £R{x))v k (16) 
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The coincidence limits Vk(v,v) up to normalization coincide with the Hadamard- 
Minkashisundram-De Witt-Seely coefficients [2J [H]. In particular, on RW space 
time, Vk((t,x),(t,y)) depends on x and y only through z 2 = (x — y) 2 and the 
coincidence limit Vk((t,x), (t,x)) of the t>& depends only on t. The same hold true 
for derivatives of Vk wrt z 2 . 

The regularization at light like separated x, y is done by adding an ie(x° — y°) to 
<t(x, y), letting e \ and taking the real part. We do not need this in the following, 
as we will approach the coincidence limit from spatial directions exclusively. 

It is shown that for n G N, the expression [{4>{x)4>(y)) w — $) n (x,y)} can be ex- 
tended to a function in C n (C v x C v ) where C v is some convex normal neighborhood 
of v. The quantum averaged field fluctuations are defined as 

<0 2 ) = Km [(^(y^-Sj^y)). (17) 

> (x,y)—t{y,v) 

where n = suffices if no derivatives of this quantity are required. For 77 7^ 0, the 
corresponding term in fllOp in the energy momentum tensor vanishes classically, but 
not quantum mechanically. For rj = 1/3, the so-defined energy-momentum tensor 
is conserved. For RW spacetimes and states of the form fl5]), the quantum averaged 
energy momentum tensor depends only on t and is diagonal so that only energy 
conservation is non trivial 

T$ 3) ) =p + H(p + p) = 0, (18) 

where p = ( * ) and p = a 2 (T- - A -* ) are energy density and pressure, 
respectively, see e.g. [19| . 
Moretti also showed 



9" V (T^Xx,, = P-3P 

1 



1 ' - 2 



+ -^rVx + cm 2 + c'm 2 R + c"DR (19) 

47T 2 

where v\ = V\(t, m 2 ) = v i((t, v), (t, v)) and c, d ', c" can be calculated from A, B, C 
and D. The D'Alembert operator here does not require point splitting and hence on 
flat RW space time can be replaced by d 2 jdt 2 + 3Hd/dt. Interestingly, this equation 
can be seen as the equation of state and the right hand side gives the deviation of 
quantum matter from the state equation of hot (relativistic) matter p = p/3. For 
Minkowski space and u the Minkowski vacuum, we expect p = p = which can be 
achieved through A 2 = Ae" 2 " 1 /m 2 with 7 the Euler constant. V\ has been calculated 
in [3] for the case of flat RW space time as 



v x = — (HH 2 + H 4 ) +— _-f OR (20) 
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In the case of a massless field with conformal coupling m? = 0, £ = 1/6, the last 
term on the right hand side of (1191) is called the conformal anomaly. 



Combining ([9]) and (TT9]) one arrives at the following equation of motion 

I \ f d 2 .. d 



-R = 8ttG 



4tt 

d 2 



1 v 1 + cm 4 + d m 2 R + c" DR } (21) 



We would like to find an expression for + 3iJ^) (0 2 ) w A ^ that does not include 



terms ( 00 } since only then second order time derivatives of the field do only 

occur on the left hand side of fll]). It is clear that j| (4 >2 ) ul \£ equals {<t>(j))u,\£ + 

(#) Wi a, € . Let fc(t) = /( s , a ')|^=t. Then, = (d 2 + d 2 s , + 2d s d s ,)f(s, s')\ s=s , =t . 
Applying this to our problem and using the fact Lemma 2.1] that 

<0(D-£ J R + m 2 )0)^ = ^ Ul , (22) 



we conclude that 



+ m Jt) {4>2)uJ ^ = + 2 M "** + ^ ~ m2 )^>^A, 5 + J^i" ( 23 ) 

Here A stands for the Laplacian on M 3 . We summarize the discussion in the following 
theorem: 

Theorem 2.1. The equation of motion for semi-classical Einstein equation on flat 
Robertson- Walker space time can be written as follows: 

-R = 87rG{(6£-l)((0 2 )^ + a- 2 (0A0) WiAC ) (24) 
+ [(2 - 6£)m 2 - (1 - 60£R] (<!>\ M + ™\~ 2 + cm 4 + d m 2 R + c" OrX 

with R = -6(H + 2H 2 ), DR=(^ + 3H%)R and v 1 given by ffijjj. 

From the above analysis it is clear that we need to compute (0 2 ) w A g and its 

second time derivative in terms of the functions G^. . . and a,H, H , R, R, R . . . in 
order to combine (jH]) and §M§ to a closed system of equations. 

Hence we need to calculate those terms in Sj(x,y)\ x o =y o =t , AgSj(x, y)\ x o =y o =t and 
dx° y)\x°=y°=t that either are singular or contribute a (time dependent) quan- 

tity to the spatial coincidence limit z— (x — y) — > 0. 
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3 Leading Terms of the Hadamard Distribution 



We start the task described at the end of the preceding section via a perturbative 
calculation of the singularities in the leading term S)'(x,y) = — (47r 2 )~ 1 of the 
Hadamard distribution. In normal coordinates X around y, a takes the form 

2a = r)^X^X v (25) 

with i]^ u the Minkowski metric. Thus the computation is reduced to finding normal 
coordinates perturbatively in dependence of the canonical coordinates. We will need 
this expansion to fifth order to fix all singularities of \im$)'(x,y) and its first and 

second derivatives as well as all homogeneous terms. The normal coordinates have 
been calculated perturbatively to fifth order in [1] (using computer algebra support), 
formula (11.12) to (11.16) for a general metric. We plug in 

= S^S^a 2 + 5?(6°6; + fyfyH™ (26) 

where denotes the nth derivative of H with respect to time and the L n are 

L = H (27) 
Li = H + 2H 2 (28) 
L 2 = H + QHH + AH 3 (29) 

L 3 =H + SHH + QH 2 + 2AHH 2 + 8H 4 (30) 
Using these formulae we get 

2a = z 2 - a 2 F 2 |l + Hz + ~(H + H 2 )z 2 Q + ^H 2 a 2 f 
+ ^{H + 2HH)4 + ^(HH + 2H 3 )a 2 z?z 
+ jgg(3ff + 6M + 2H 2 - 8HH 2 - AH^z* 

+—(9HH + 8H 2 + 7AHH 2 + AZH^cfz^zl 
360 V ^ o 

+^(3M S + 4# 4 )aV) 2 } + 0(z 7 ) (31) 

where all time dependent terms are evaluated at y and we use the abbreviation 
z = x — y. 

Lemma 3.1. a is symmetric under exchange of x and y 

Proof. The proof is done by straightforward calculation, however on has to keep in 
mind that a, H and its derivative all have a suppressed argument yo such that they 
have to be Taylor expanded. The terms of even orders in z do not change signs so it 



7 



has to be shown, that they are not affected by the terms from the Taylor expansion 
of the lower order coefficients. Using abbreviations a y = a(y ) and the likes we get 
to the relevant orders 

4=1 + 2H y z + (H y + 2H 2 y )z 2 + \(H y + QH y H y + AH 3 y )z 3 
a y 6 

+ 1(4 + 8H y H y + QH 2 y + 2AH y H 2 y + 8H±)z* (32) 
-hAzo = - H y z - (H y + 2H 2 y )z 2 -\{H y + QH y H y + AH 3 y )z 3 



V 



l -(H y + 8H y H y + 6 H 2 y + 2AH y H 2 y + SH^ ( 33) 



-{H x + H 2 )^z 2 =-(H y + H 2 y )z 2 + -(H y + AHyHy + 2H 3 y )z 3 



y 



1 

( 

+ \{H y + QH y H y + ml + UHyH 2 y + AH 4 y ) z A (34) 



1 •■ -a 2 1 ■■ 

~T^( Hx + 2H x Hx)~^z = - —(H y + 2H y Hy)z 



1(4 + AHyHy + 2flJ + AH y Hl)zl (35) 



= I2^ + l^ yHy + 2 4) z ° 



+ 1(4//, + 4 + lO/j,// 2 + ml)zl (36) 
~(H X H X + 2H*)^z = - 1(4^ + 24)^o 



12 

T 3\ u x , 1 



- l(ff s ff w + 4 + WHyHl + ml)zl (37) 
Adding up the first four and the last two terms one gets 

1 + H y z + \(H + H 2 )z 2 + l(tf + 2HH)z 3 (38) 

±Hl + ^(H y Hy + 2Hl)z (39) 
which proves the claim. □ 

To calculate u we use the first equation of well know Hadamard recursion which 
is deduced from the Klein Gordon equation 

2g^{x) {d?<j) (<9» + (□ V - A)u = (40) 

where one has to take into account, that the derivatives are with respect to x and 
the metric is evaluated at x in the above formula, which leads to additional terms 
in the calculation. We make the ansatz 

u = 1 + pz 2 + vz 2 + pz% + rz 2 z + 04 + *l>?4 + X(^) 2 (41) 
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where all coefficients are evaluated at y . The first term being 1 and the absence of a 
first order term follow from the requirements that u(x, x) — 1 and u(x, y) = u(y, x). 

The calculation of the coefficients then yields 

u = l-±(H + H 2 )zl + ~^(H + 3H 2 )a 2 z* 

-l(H + 2HH)4 + ^r(H + SHH + 6# 3 )a 2 i% 
o 24 

+—(-18H - 36HH - 17H 2 + 38HH 2 + 19H 4 )zi 
480 

+ ^-(3H + 26HH + 17H 2 + 52HH 2 + lH 4 )a 2 fz 2 
240 u 

+^(4M + 3i7 2 + 36M 2 + 29H 4 )a\z 2 ) 2 + 0(z 5 ) (42) 

which can be checked to be symmetric under interchange of x and y by a straight 
forward calculation. This equation has been derived using computer algebra support, 
but has been validated up to second order by hand calculations. 

From this we obtain 

Theorem 3.1. The most singular order of the Hadamard distribution is given by 



- -^—(12HH + 9H 2 + 86HH 3 + 51H 4 )a 2 f + 0((z 2 ) 2 ) 
1440 

(43) 

-4n 2 W(z)\ z o =0 =JL - Uh + AHH) + 0(f) (44) 
a 2 9 d t ,, t m 2 H 2 



dx° dy° ^(z 2 ) 2 a 2 ? 2 



1 (4H + 16HH + 27H 2 + 30M 2 + 17H 4 ) + 0(&) 



240 

(45) 

The proof is done by straightforward calculations where the occurring fractions 
in powers of Z = z 2 are expanded in a power series like 

L + MZ + NZ 2 L MP-LQ^ NP 2 - LPS - MPQ + LQ 2 ~ lA . 
P + QZ + SZ* =P + ^^ Z + pa Z ( 46 ) 

Even if placeholders are inserted for the coefficients of the powers of Z and z in 
a and u the calculations remain very lengthy, as the coefficient's derivatives have to 
be taken into account. Therefore we will not show them here explicitly. 

We calculate the sub leading terms of the Hadamard parametrix in position space 
to the order relevant to calculate the homogeneous term of the energy momentum 
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tensor. This means we use the Hadamard recursion 



2g^(x)(d*a)(dZv Q ) + (DV - 2)v = -(0 X + m 2 - £R(x))u (47) 
2g'"{x){d£a){dZv 1 ) + (D x a) Vl = -(□* + m 2 - (R(x))v (48) 



to calculate 



and 



v ° = ~KG ~ v + y + ^ + AkH)z " 

+ 2i0 (J 21 ~ 120 ^ + ( 87 ~ 480 O™ + (54 - 300O^ 2 

-(76 - 540OM 2 - (58 - 360O# 4 + 30m 2 (H + H 2 ])z% 

+ 2i0 ( ~^ + ^ ~~ m &HH + (6 - mi)H 2 + (76 - 540£)M 2 

+ (58 - 360O^ 4 - 10m 2 (H + 3# 2 ))a 2 ^ + C(^ 3 ) (49) 

«i = i^o (C 1 - 50^ + ^(1 " 60 2 ^ 2 " 2{H + H 2 )H 2 + 5(1 - e^Rm 2 

+ 15m 4 ) + 0{z) (50) 



4 Mode Expansion 

Two things remain to be done: The singular terms in f)13p have to be related to the 
mode expansion etc. of the two point functions and the homogeneous term of 

the energy momentum tensor has to be calculated. 

First, the expressions in Theorem 13.11 can easily be Fourier transformed with a 
method explained e.g. in [2j. Using the distributional Fourier transforms 
^(\z\~ 2 )(k) = 2ir 2 k~ 1 , JP(\z\~ A )(k) = —ix 2 k with k = \k\ and ignoring for now 
homogeneous terms that will only contribute terms proportional to the delta distri- 
bution in the zero mode and positive order terms in Z, we get the mode expressions 

^($)'\ z0=0 )(k) = a- 2 ^ (51) 

# (jy|*o =0 ) (k) = -a- 2 ^ (52) 

^ ( d d „ \ /T * N A k H 2 , . 

*\WW^^) {k) = a 2 +a ^ (53) 

Here it needs to be taken into account that all singular orders in z 2 that are 
encountered are locally integrable functions on M? with the exception of the \ term 
in d^a^'l^=o- Continuation of the ^s^s[{<j>{x)<f>{y)) - Sj n {x, y)]\ z o =0 for n > 2 
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to z = thus implicitly induces a rotational invariant regularization of -4, as it 
can be written as a difference of already regularized terms and the newly found 
continuous function. The analytic regularization [2J of \z\^ is well defined except 
for poles at C — ~~ 3, —5, —7, .... It is clearly rotational invariant. Furthermore \z\^ 
with £ = —4 is the only invariant extension of on M 3 that preserves the scaling 
degree, as there are no rotation invariant linear combinations of first order derivative 
of the delta distribution at zero. But the difference of two extensions with the given 
properties has to be such a linear combination, [2J. Hence we conclude that it has to 
coincide with the regularization induced by the continuation of -£j}-^o[(<P( x )4 > (y)) ~ 
S) n (x,y)]\ z o =0 to z — 0. The Fourier transform (153]) can thus be calculated based on 
the analytic regularization prescription [2j. 

We now turn to the mode expansion of the sub leading terms. Expanding the 
logarithmic factor with the help of ( |3T|) yields 



log 



<r(x,y) 
A 2 



oc 



n=l 



log(^J+log(^)-^(--i/ 2 aV + 0(^)) , (54) 



where the infinite sum can be truncated after a few relevant terms. Inserting (131)) 
to expand the powers a(x,y) k in (fT3|) . we easily see that the singular contribution 
from the second term in (fT3"j) to [(<p(x)<p(y)) — S) n (x,y)]\ z o =0 is 



1 „. _ ,. _ (, ( a 2 \ , , „, 



v ((t, x), (t, y)) log — + log(^) . (55) 



4tt 2 

Likewise, using ^U^o =0 = Y^=o(h,H 2 a 2 z 2 + (^(z 4 ))™, it is easily shown that 
the expansion of derivatives of the second term in (|T3|) with respect to and d ^t dy0 
at equal time into singular orders in z 1 which is truncated after some sufficiently 
high order contains only terms ~ log(z 2 ) and ~ (^) fc , k > 0, where the coefficients 
are made of log(A), log(a) and Uk((t,x),(t,y))-terms and their derivatives w.r.t. 
z 2 along with their time derivatives; for the case with the two time derivatives an 
additional term with the singularity structure \ needs to be taken into account as 
well. 

One could calculate the coefficients of the singular terms described above by 
calculating the derivatives of (T4*9|) and (15 0|) and Fourier transforming the result by 
hand. However these calculations are very tedious, therefore we proceed differ- 
ently here. As we have analyzed the general shape of the terms to account for, 
we can Fourier transform these terms and include them into the mode expansion 
with up to now undetermined coefficients. Doing so, we use that jF(log(|2:|))(fc) = 
— 47r 3 / 2 r(|)A; -3 . Note that this distribution at does not extend a locally integrable 
function and hence requires a regularization prescription. The details can be found 
in appendix [A] The remaining singular orders that occur have mode expansion 
^(^log(|^|))(A?) = 24vr 3 / 2 r(f)Ar 5 and ^((r^)")(fc) = (2n) 3 (-A) n 5(k), where we 
again neglect the latter for the time being. 

Next we multiply the leading terms in the singular order expansion with the 
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powers of a to get expressions similar to the two-point function 



a~ 2 

= -^f ( 57 ) 
. a 2 k a^H 2 . . 

* - = - + ~w (58) 

So, taking now into account the discussion above, we consider the ansatz 

= + ^ + ^ + (60) 

that captures all possible singular orders except for those concentrated in the zero 
mode. 

Note that S) n fulfills the Klein-Gordon equation up to terms that vanish in the 
coincidence limit and a zero mode term. Plugging our ansatz of the singular order 
expansion into the system of equations (jH]) we obtain a couple of equations that help 
us to determine the unspecified coefficients 



A3 As _ 2a~ 3 /3 3 2a~ 3 /3 5 

2k 3 2k 5 ~ 2k 3 2k 5 [ ' 

a(H + 2H 2 ) /3 3 (-£H + ra 2 )a + aa 3 + m 2 )a 3 a 3 + aa 5 



2k 2k 3 2k 2k 3 

+ -m (63) 

2a 4 #(# + 2H 2 ) + 71 73 -2(-gfl + m 2 )a 4 H + 2a/3 3 
2/c 2/c 2/c*^ 2/c 

2(-£ J R + m 2 )a 3 /3 3 + 2a/3 5 
2P 

where we suppress undetermined orders. 

Assuming these equations to hold order by order in k we get two sets of equations 

a 3 = 2a~ 3 [3 3 (65) 



(64) 



71 = -2a/3 3 + 2a 4 H ( ( - - £j + m 2 J (66) 
«3 = «" 4 7i-Q-e) J R-^ 2 (67) 
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and analogously 

a 5 = 2a _3 /3 5 (68) 

73 = -2a/3 5 - 2a 3 /3 3 (-£R + m 2 ) (69) 

a 5 = a~ 4 7i — a _1 /3 3 — a 2 a 3 {—£R + m 2 ) (70) 

We now want to solve these differential equations. The equation ( I65p + (l6"6"j) -& (l6"7j) 
leads to ^ 

7l - 2Hj = j Q - £^ (i? + 2HR) + a 4 #m 2 (71) 

which can be solved using standard methods. The solutions for and /3 3 can then 
be calculated straightforwardly without solving differential equations. 

Lemma 4.1. The solution to the system / fg3]) -/ fg7|) is 

i]R + m 2 ) + Aa~ 2 (72) 

(73) 

f Aa 2 (74) 







«3 = 










a 3 




~ 4 










7i = 





where A is some constant. 

Proof. A differential equation of the form 

x(t) - 2/(*)x(*) = g(t) (75) 

has the solution 

x(t) = J g(t') exp I 2 J f{t")dt" \dt' + A exp I 2 ^ /(t')dt' (76) 



o 



*2 

with A being constant. In our case /(*) = H, such that exp I 2 / H(t')df 



We get the solution 

f 



a 2 



o 

i 

1 



--i)R + m 2 )a 2 



o 



a 2 (fl) 
a 2 (t 2 )' 



li = a 2 j Q Q - £J (i? + 2F^ + #m 2 ^) + Aa 2 (77) 



dt + Aa 2 (78) 



y (Q-e)^ + m 2 )+Aa 2 (79) 
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where the constant A absorbs all constant terms and thus changes from line to line. 



2a/3 3 = -71 + 2a A H ( Q - Cj R + m 2 J (80) 

= -2a' H A R + m ^ - ^ (l- ^ R - 2Aa 2 H (81) 

+2a 4 #(Q-^) R + m 2 ^j (82) 

" ' / ' Z)R- 2Aa 2 H (83) 



2 V6 



and 



« 3 = «" 4 7i - Q - c) ^ - rn 2 (84) 

= K Q " * + m ') + -(l-t) R - m2 (85) 

=_ KG~ e ) /?+m2 ) +Aa " 2 (86) 

yield the other two coefficients without solving differential equations. □ 

Applying the same strategy as before we get the solutions for the second system 
of equations 

« 5 = y| Q-e) (R + 5HR-HR) -3^-A^R 2 (87) 



- (4if + QH 2 + 6£i?) m 2 - 3m 4 



a 3 

As = y"5 (89) 

73 = ~| Q - e) + ^ + ™) - Q - e) ^ 2 (90) 

-2(/7 2 + ^)m 2 -m 4 | (91) 

where we suppress the integration constant as it vanishes for a pure state due an 
argument given below. 

Theorem 4.1. The singular part of the Hadamard modes is equal to the singularities 
of the two-point function of a pure state if the integration constant A in Lemma \Ji~l\ 
vanishes. That is, 
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• (%) Vk : fi^kfi^k - $? (<j> ^ k = | + O (k~ 6 ) 

• (u)Vk: Sj^k + tf^fc > + O (k~ s ) 

holds if and only if A = 0. 

Proof. The first claim is equivalent to 

a 2 H 2 + a~ 2 7i + a 2 a 3 - a 2 H 2 = (92) 
gT 2 7 3 + a 2 a 5 + a A H 2 a 3 + a 3 7i + 2aHf3 3 = (93) 

where equation fl92|) consists of the terms of order kr 2 and equation fl93|) represents 
the order k~ 4 . Equation (I9"2"j) gives = 0. Analogously equation (I9"3"|) forces the 
other integration constant to vanish. 

To show that the second claim is satisfied we use the fact that is a solution of 
the Klein-Gordon equation up to O {k~ 5 ). Then we use a form of the well known 
deformation argument of Wald, which is explained in [16]. We deform the spacetime 
for early times such that a is constant at early times. Thus the claim holds exactly 
at early times, as we just get the Minkowski two-point function. Then we propagate 
with Equation ([6]) to some later time, where the spacetime is undeformed. Thus 
the claim is still satisfied at that time because the first claim shown above ensures 
that Sj^kfj^tk > + O {k~ 6 ). Finally, we undeform the spacetime at early times 
and can thus show the claim to hold for all times using Equation (jSJ). □ 

From Theorem 14. II we obtain a (tentative) parameterization of initial conditions 
of the field degrees of freedom for the equation of motion for pure states. For k > 1 
we set G^k = fi^k + a (k), G^ tk = $)^),k + b(k) and G^ jk = {\ + G^^/G^k 
where a(k)k 7 and b(k)k 7 are C£°((l,oo)). Furthermore, a(k) needs to be choosen 
such that (j00,fc > for all such k. This can always be achieved as the leading term 
aT 2 /k is positive. For k < 1, these functions need to be continued to C£°([0,oo)) 
such that the relations G^^ + G W7Ti k > and (JSJ) are preserved. Obviously with 
such initial conditions the trace of the initial energy momentum tensor is finite. In 
particular, Hadamard states fall into this class. 

From the construction of Sj^^, $)(<j>n),k an d ^^^^ it follows that these quantities 
fulfill ([6]) up to order k~ 7 and k~ 5 , respectively. Thus, if the differece of and 
Sjt/t^k (multiplied by k 2 ) is integrable for large k, this property will prevail after 
an infinitisimal time step. The argument for the difference of G n7T ^k arid S)^^ is 
analoguous. 

For non pure states, a positive c(k) can be added to G n7T k such that k 5 c(k) G 
C fe ([0,oo)]. 



15 



5 Equation of motion 



We begin with the easier case where the field is conformally coupled to the mean 
curvature, i.e. £ = 1/6. In this case, ( I24p simplifies considerably, as it only contains 
the field dependent term (00) Wi a,£- This term, up to zero mode contributions, is 



^ 3 Jm.3 



a 2 m 2 



G<t>4>,k — — — \- 



2k 4k 3 



dk, 



(94) 



since a 3 = — The renormalization prescription of J R3 k 3 dk at is given in 
Appendix [A] 



The zero mode terms in the conformally coupled case are 



1 1 (a 2 

y 2 R- log i Y2 



(95) 



where the first term stems from the first equation in Theorem 13.11 and the second 
one from (155]) . In the conformally coupled case, Vq(x, x) takes the simple form 



v = — 



m 



(96) 



Wrapping up, we obtain 

Theorem 5.1. For the case of conformal coupling, the equations of motion is 



-R 



87TG 



m 



<,7T J 



a 2 m 2 
~~2k + 4k 3 



,ik - nr ( + c \ n 



240 



-m 



dn 2 \ J V 288C 



n , , c" )UR + — log — 

7T 2 / 47T 2 V A 2 



32 7T 2 



where the k 3 integral is regularized at zero as explained in AppendixlM 



(97) 



Wald's fifth axiom [XT] [T8] states that the equation of motion should be a second 
order equation. (19T|) is second order if and only if the renormalization constants 
A,B, C and D in ffTT]) are chosen such that c" = — l/(28807r 2 ). Note however, that 
for this value of c" the joint system differential equations ([6]) and (l97j) is an implicit 
(infinite dimensional) system of differential equations and for c" ^ — 1/ (28807T 2 ) it 
is explicit, which is of advantage for proving existence and numerical solutions. 

The log(a) term that seems to be missing in [13] is due to the fact that the 
infra red prescription we employ for J R3 ^grdk (cf. Appendix A) is not a dependent, 
whereas the prescription used in [13j is, leading to the absence of a stand alone 
log(a) term in the latter case. 
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For £ 7^ 1/6 the homogeneous terms of the Hadamard parametrix take a much 
more complicated form, as second derivatives of the parametrix must be taken into 
account. Ignoring a prefactor |^ we have to calculate zero mode terms of 

£t,m =(6£ - 1) {d Xo d yo + a' 2 A) (-2^ - v log (2a) - ± Vl 2a log (2a^j 

- 1 [(2 - 6£)m 2 + (6£ - l)£i?] (if + 2tf 2 ) 

+ log (^j ^(6£ - 1) (9 X0 9 W + a~ 2 A) - ^i2^ 

- [(2 - 6£)m 2 + (6£ - t, j (98) 

in the coincidence limit x — > x'. This can be done using the expansions of u, Vq, 
v\ and a 2 given in Section 3, applying the differential operator, setting z° = and 
extracting terms of z-order zero. With a computer aided calculation this yields the 
result 

= - — (ffl + 53HH + 11H 2 + U1HH 2 + 3H 4 ) 
30 V / 

+ ^ (9H + 113HH - AAH 2 + 11HH 2 - 277 H 4 ^ 
- 6£ 2 (h + 12HH - UH 2 - 38HH 2 - 68# 4 ) - 108£ 3 (h + 2# 2 ) 2 
((21i? + 44# 2 ) - 6£ (3H + 8# 2 ) - 36£ 2 (h + 2# 2 )) 



2 

m 

+ T 

4 

m 



+ log (J^J + 21HH + 42if 2 + 152M 2 + 116# 4 ) 

+ ^(4H + 28HH + 31H 2 + 106HH 2 + 58# 4 ) 

- £ 2 ^3ff + 21M - 6H 2 - 36HH 2 - 72# 4 ) - 108£ 3 (h + 2# 2 ) 2 

+ \m 2 (1 - 6£) (H + 2H 2 ) - m 4 (1 - 3£) ] (99) 



Although this result looks very tedious, one can simplify the result subsuming 
terms proportional to DR, and m 2 R into the renormalization degrees of freedom 
and the H independent term into a redefinition of the scale parameter A. 
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The simplified form is 



Z im = -— \--UR - —R - —R 2 + 13HH 2 + 23H 4 
^ 30 V 3 6 36 

+ ^ (--OR - | R - ^R 2 + 23HH 2 + 43H^j 

(--OR --R- -R 2 + 2HH 2 + 4H 4 ] - 3£ 3 R 
V 6 6 2 J s 



+ ^ ( {~R + 2H 2 ^j - 6£ (-±R + 2H 2 ^ + 6^) + ^ (1 - 60 

+ log f^V _ 1 f-Ini? + ^ - 4M 2 - 4tf 4 
S \\ 2 ) I 60 V 2 6 

5^ V 3 12 

- £ 2 ( ~DR - Ir 2 ) - 3£ 3 R 2 + -m 2 (1 -6£)R- m 4 (1 - 30 I (100) 



Plugging in the calculated homogeneous term we obtain the equation of motion 
for the general case. 

Theorem 5.2. The equation of motion in the general case is given by the following 
expression 



- R = 8nG U6£ - 1 



G-n 



^ 3 JR3 



k 2 G 



<j><j),k 



7T1T,k 

a~' 2 k «3 «5 
~~ 2 2k ~ 2P 



a 2 k a 4 H 2 + 71 73 
2F 



dk 



dk 



G 6 : 



a 0:3 
~2k ~ 2k? 



dk 



I 35£ _ 5 

-Xt.m + 7 n ^1 + cm 4 + cWi? + e"OR 



47r 2 ^' m ' 4tt 2 

TTie regularization of the integrals over k~ 3 at zero is given in Appendix L"4l 



(101) 



We would like to point out that due to the 4th time derivatives in 73 and a^, (11011) 
again leads to an implicit system of differential equations, as the leading order time 
derivative can not be isolated, regardless of the renormalization degrees of freedom 
and the specific form of (E^ m (t). Also, Wald's fifth axiom can never be fulfilled in 
this case. 
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A Regularization prescription on k~ 



Here we give the details of the regularization of the seemingly infra red divergent 
integrals 



a 2 m 2 
" ~2k + 



(102) 



and related. For any Schwartz test function % with x(0) = 1 This integral can be 
re-written as 



m 



G(t>4>,k ~ 1" TTo 



2k Ak 3 



:i-x(*))+ 



G<t><t>,k — 



2k 



m 



X (k)dk+ I -^X(k)dk (103) 



Here the first integral is regular. The second - seemingly infra red divergent - 
integral needs to be properly traced back to its definition as distributional Fourier 
transform of log(r) smeared with the test function x, which gives a regular integral 
after shifting the Fourier transform to the test function. 

Let us quickly calculate the Fourier transform in the sense of tempered distri- 
butions of the locally integrable function log(r), r = |x| for the convenience of the 
reader. Let <p(x) = ip(—x), then 



(^(log(r))(A;)^> 



(log(r),^)) = hm^<^,^)> 



lim — 

C/o d( 



2 C+3 7r 3/2 



em 

r(4) 



(AT*- 8 ,*,) 



(104) 



Since £ approaches from below, fc -l » -3 is locally integrable. Furthermore, (k~^~ 3 , cpS 
can be analytically continued in ( to C\No [2]. The Laurent series of k~^~ 3 in £ = 
has a pole proportional to 5 . Hence, for (p(0) = 0, (fc - ^ -3 , (p) is analytic at ( = 
and in particular has a finite £ derivative. So let us from now on assume (p(0) = 
and perform the ( derivative at the right hand side of fj!04[) . Due to the pole of 
r (— |) at £ = 0, the only term that contributes in the limit £ /* is 



lim2^ 3 -V/ : 



: iW)r'K) 
r(-|) 2 



<^- 3 ^> = -4 7 r 3/2 r(^)<A ; - 3 ,^>, (105) 



where we used the Laurent series expansion of T at zero. Since, by dominated 
convergence, (k~ 3 ,(p) = J R3 k~ 3 f(x)dx, we see that J^(log(r))(/c) is a regularization 
of -47r 3/2 r (|) k~ 3 . Let now ip(0) ^ 0, and x(k) = e~ fc2 . Since <p{x) = [ip(x) - 
</9(0)e~ fc l 2 \ + (p(0)e~ k I 2 = tpi(x) + (p(0)e~ k I 2 . We have already shown how to 
evaluate J^~(log(r))(/c) on ip\. It remains to calculate 



^(log(r)),^(0)e- fc2 / 2 



y(Q) 

(2vr) 3 / 2 



log(r),e 



-k 2 /2 



y(Q) pi 

(2tt) 1 /2 ^2 



(106) 
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